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The consequence of jet interacting with a warped disc 
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ABSTRACT 

The sprinkler pivots on a bearing on top of its threaded attachment nut. It is driven 
in a circular motion by a spring-loaded arm pushed back by the water stream which 
returns to "impact" the stream. The water stream can thus rotate around a fix axis. 
Analogously in our universe, the outflow or jet formed by the relativistic plasma 
corresponds to the water stream of a sprinkler; and the baryons in the tilted outer 
accretion disc or torus play the role of "impact arm" . Then the jet aligning with 
the inner parts of a warped disc can directly "touch" the outer region of the disc. 
The resultant collision between such rapid leptons and slow baryons automatically 
accounts for the main features of broad-line region of active galactic nuclei. Moreover, 
it naturally provides a channel of dissipating the angular momentum of an accretion 
disc, which has long been a problem in theory of accretion disc. 

Key virords: accretion, accretion discs — (galaxies:) quasars: general — 

galaxies:individual (3C120,3C279,3C273,3C345,OJ287) 
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1 INTRODUCTION 

An active galactic nucleus (AGN) is comprised of three ma- 
jor elements, an accretion disc, a central spinning black hole 
(BH), and a jet. The complex interactions among them re- 
sult in fascinating phenomena of all bands, from radio, op- 
tical to X-ray and Gamma-ray, and displaying in all forms, 
from images to spectra and time varying effects. Although 
tremendous progresses have been made in understanding the 
mechanisms underlying them, many fundamental questions 
are still open. 

Firstly, frame dragging produced by a rotating BH with 
angular momentum S causes precession of a particle if its 
orbital plane is inclined in relation to the equatorial plane 
of the rotation object. This is known as the Lense-Thirring 
effect l|Lense fc Thirrin3ll918h . with a precession velocity, 
flLT oc of strong function of radius, r, so that the inner 

parts of the disc tends to align with the direction of the BH 
spin. 

In the case where the BH spin is misahgned with the 
rotation direction of the disc, the disc becomes twisted 
and warped. Such a warp can propa gate in the d isc in 
a diffusive manner through viscosity (|Pring:lel Il992h . And 
the warp can make a l arge increase in the accretion rate 
l|Lodato fc Pringlell2006l h 

On the other hand, when a BH is in a binary sys- 
tem, the inner regions of the disc remain aligned with the 
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spin of the central BH. Whereas, the outer region pre- 
cesses around the orbital axis, owing to the tidal torque 
of the com panion object (|Katz 19731: Sillanpaa et al.l 1988 : 
iKatzil 19971 : lRomero et~al]|2000l : lcaproni fc Abrahamll20o'4b : 
iMartin et al.l 12009^ 1. This results in a warped and twisted 
disc. 

The apparent lack of correlation between the orienta- 
tion of the radio jets and the plane of the host galaxies 
disc m ight also be attributed to warped discs (jSchmitt et al.l 
I2OO2I) . 

Then there is a question: what if a disc is so warped 
that the jet axis aligning with the inner disc touches the 
outer disc? Once such a warp induced jet-disc interaction 
occurs, the inner and outer parts of a warped disc interact 
in a new way other than the viscosity, which influences the 
precession velocity of jet, accretion rate and stability of the 
disc. 

Secondly, accretion disc has been proven to be one effec- 
tive mechanism of extracting gravitational potential energy 
and converting it to radiation powering AGNs. No matter 
a disc is warped or not, the accreted material at the inner 
region of the disc has to dissipate its angular momentum in 
order to fall on to the BH. If there is no external torque, the 
angular momentum loss is achieved through the accretion 
disc by outwardly transferring angular momentum. Accre- 
tion onto an AGN black hole will be limited largely by how 
rapidly angular momentum can be lost, which determines 
the power of AGN directly. The angular momentum trans- 
ferred to the periphery must have a reasonable way of loss. 
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or the accumulation of external angular momentum would 
cause instability in the accretion disc. 

Thirdly, one of the fundamental and defining charac- 
teristics of AGNs is the presence of broad emission lines in 
their optical and ultraviolet (UV) spectra. The investigation 
of the profile and intensities of these lines show that they 
originate in dense gas moving very fast around the central 
engine, which are thought to be fully exposed to the intense 
radiation from the engine. 

The outer region of disc and wind from disc as 
possible origins of B LR have been studied extensively 
l|Eracleous et al.ll2009l ). However, the complexity of obser- 
vations, e.g., the double-peaked broad emission lines are 
thought to originate in disc, while the local viscosity heat- 
ing is usually insufficient for the observed luminosity of the 
emission site; and the correlation of BLR (optical and UV) 
with Fe Ka (X-ray) lines indicate that the origin of BLR is 
still an open question. 

Interestingly, each of the three questions, the preces- 
sion property, transfer of angular momentum, and origin of 
BLR, appears to be unrelated and difficult, but all of them 
together point to a simple scenario: sprinkler, with the wa- 
ter stream replaced by the jet and "impact arm" replaced 
by the outer parts of accretion disc. 

The collision between jet (formed by relativistic elec- 
trons) and the outer disc (slow baryons), automatically 
kicked out baryons from the disc, which provides a chan- 
nel of angular momentum dissipation of a disc. Moreover, 
such a jet disc interaction influences the overall precession 
velocity of the disc. Finally, the baryons in the disc, carrying 
velocity gained through the collision with electrons, spread 
around the collision site, which naturally account for the 
major features of BLR. 

Therefore, sprinklers in our yards, may provide a clue 
to understand what happens to "sprinklers" of the garden 
of universe. 
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The accretion rate of mass M can be expressed by, 

-S = f [l-(f )^/^ 

where R^, — 2G Mp / is the Schwarzschild radius of a BH. 

Substitute the Kepler angular velocity, Q, — ^k{R) — 
{GMp/R^y^^, and Eql2]into Eq[T] the angular momentum 
loss rate can be written in terms of M, 



i~{-2.7?^-i[l-(^)^/^]-^(^^)^/^}| 
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(3) 



Such a L is the angular momentum required to be trans- 
ferred to outer disc per unit time in order to maintain per- 
sistence power of an AGN fueled by M . There must be an 
efficient channel of angular momentum dissipation, other- 
wise the piled up angular momentum in the outer disc can 
cause instability of the disc and stops the accretion process. 

The sprinkler scenario provides a simple answer. The 
collision between jet and disc results in scattering matters off 
the disc, which should be in Kepler motion around BH and 
hence carrying angular momentum. Such a mass loss rate of 
the outer disc can be estimated by, rhioas = MdiscL/ Ldiac, 
where Mdisc ~ nrup and Ldisc ~ nLp are mass and angular 
momentum of the disc respectively (n is the total number 
of baryons in the disc) . 

If we use the angular momentum of the outermost ring 
to represent that of the whole disc, the mass loss rate reduces 
to 



mpL ■ . ,Rt -.1/2-1 
Lp ro 



(4) 



where Lp = n^p ujq = mp(G Alp tq )^^^ is the angular mo- 
mentum of a proton. Thus, a channel of dissipating angular 
momentum is set up, at a warped disc as shown in Fig.l. 



2 THE NEW MODEL 

Once a disc is warped as shown in Fig.l, the jet-disc inter- 
action occurs. The collision sites separate from the central 
engine by a distance ro. The opening angle of the jet cone 
is denoted, do, and the misalignment angle between the jet 
axis and the axis of the precession cone is uq, as shown in 
Fig.l. 



3 TRANSFER OF ANGULAR MOMENTUM 

At different radii of a disc, gaseous material is orbiting with 
differing angular velocities, Q, determined by Kepler's third 
law. Because of the ever present chaotic thermal motions 
of fluid molecules or turbulent motions of fluid elements, 
viscous stresses are generated. This type of transport process 
is known as shear viscosity. 

The torque exerted by the outer ring on the inner 
one (= —the torque of the inner on the outer) is G{R) = 
27ri?z/E R^ where u is the kinematic viscosity, E is 

the surface density of the disc, and Q{R) is the angular ve- 
locity. Such a torque corresponds to a transfer of angular 
momentum (per unit time) to the outermost ring. 



4 THE COLLISION PROCESS 

The jet-disc interaction dissipating the angular momentum 
can be calculated simply by the collision of relativistic lep- 
tons in the jet with the slow baryons in the outer disc, as 
shown in Fig. 2. Major parameters in the study of such a 
jet-disc collision are: the volume of the interaction zone, Vi', 
the number density of baryons (dominantly protons) in the 
outermost disc, n^; the number density of leptons (relativis- 
tic electrons) in the jet, n^, which travel at speed, v^. The 
number of leptons scattered into a solid angle between Q 
and Q dQ per unit time is given, 

dNi — 2-Ka{Q) sin QdQ ■ n^v^- Up Vi ■ (5) 

The number of leptons ejected out from the jet per unit 
time is estimated, A^e ~ S,M / nip, where ^ is a dimensionless 
coefficient, ^ = 1 means accreted matter is totally ejected 
out. With an opening angle of jet, ^ shown in Fig.l, we 
have UeVe = Ne/[Tv{^ro9o)^]- The change of linear momen- 
tum (per unit time) of leptons involving in collision is given, 
dPe = 2\/2me iJsin ^d Ni, where E = 'y^nieC^ is the en- 
ergy of a relativistic electron, with mass and Lorentz fac- 
tor 7e . For simplicity, the Rutherford scattering cross section 
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l|Goldstein et al.ll2005h . a{Q) = i(^2f^ )^ csc^ |, is used, 
where E is the energy of the scattering particle. Ahhough 
it is derived in the case of a repulsive scattering process, it 
is still a good approximation in the collision of the jet and 
disc. 

Since the collision occurs predominantly between elec- 
trons and protons, the interacting charges satisfy, (ZZ' )^ = 
1. Substituting EqO into the linear momentum change 
yields. 



dPe = V2me E{ 



2E 



.2 SCMnpV; cosf 
rriprlel sin^ | ' 



(6) 



The mean force corresponding to the interaction is 

F= dp,. (7) 

The number density of protons in the outermost ring is Up, 
each proton occupies a volume 1/ Jip , the functional rela- 
tionship between the impact p arameter and the scattering 
angle is (jColdstein et al.ll2005l ). s = ZZ'e^ cot f /(2S). In 

— 1/3 

such case, we have s ^ rip , the lower limit of the integral 
can be obtained by, 

7' „2 



-1/3 



ZZ' 



eo 



„ ■ cot 
2E 2 



(8) 



correspond to Ea llll and the solid vertical ones correspond 
to EqOin which L is estimated under the binary model, with 
radius ro eq uivalent to the outer disc radius rg under th e 
tidal torque (jPapaloizou fc Terauem|[T995l : lLarwoodlll997l ). 

The shadowed area of Fig. 3 corresponds to an accept- 
able parameters of collision process, the predicted torque of 
which can precess the jet at a peri od of tens of years, which 
is expected by Newtonian torque (lRiegeij|200i ). 

The volume of the collision region is estimated, Vi ~ 
2ro6o ■ TvR^i^g, the magnitude of which can be inferred from 
the cross of vertical lines (correspond to different scenarios) 
with the inclined lines of Fig.3. 

The shadowed area corresponds to the case that the 
lOyr precession period is driven by the Newtonian torque, 
so that the size of the outermost ring of the disc can be 
estimated, Rring ~ [Vi/ (2nro9o)]^^'^ , which is ~ 3.5 x lO^^m 
for 3C 120, and ~ 7.4 x lO^^^m for 3C 345. 

In contrast, if the ~10yr precession period originates in 
the Lense-Thirring precession in a single BH, the jet-disc 
collision results in a twisted outer ring of the disc. Because 
the collision makes the liner momentum of baryons in the 
outer disc vary as Fig. 2. Such a twisting of outer ring appar- 
ently changes with the jet precession, so that the whole disc, 
inner and outer, varies at the precession rate dominated by 
the inner disc. 



Integrating Eq[71 we get 



F = 




2E 



— ? - lW2m,E(— )^ l^CMrtpy^ ^ ORIGINATED IN COLLISION 

1/3 ^ " ^2E ^ m„rlel ' ^ ' 



Eq[9] corresponds to a torque. A'' = Fro, which makes 
the jet to precess around the core at angular velocity, Hp, 
determined by, N — Qp x S, so that, 



A'' = fJpS'sinao — 



27r(l -I- 2)5 sin qq 



(10) 



where P is the precession period measured in the observer's 
reference frame, and S = GMp\at\/c is the spin angular 
momentum of the primary BH, with a, its dimensionless 
angular momentum. Such a precession resembles a sprinkler, 
with water stream and impact arm replaced by the jet and 
the baryons at the outer disc respectively. 

Thus, according to the collision process described by 
Eq[5] to Ea llOl the mass loss rate, M, required to precess 
the jet at the period, P, is, 



As a result of the jet-disc collision, baryons in the outer disc 
scattered by high-speed leptons in jet spread around the 
collision site. 

The momentum and energy exchange between slow 
baryons and relativistic electrons are shown in Fig. 2. The 
velocity distribution of scattered baryons can be calculated 
simply by the conservation of momentum. 



v = 2\ 



-ccosS, 



(12) 



M 



7r(l -f z)mpro9oGMp\a,. 



sm Qo 



(11) 



1] 



where p — /{2E). Fore simplicity, ^ = 1 is assumed in 
the calculation of Ea llll Notice that M of Eg 1 1 1 1 represents 
mass loss rate of the relativistic jet, which is thought to be 
comparable to the accretion rate of Eq[3] derived from the 
angular momentum transfer of a disc. 

The parameter spaces of M , ro and Vi corresponding to 
the two mechanisms given by Eq[3] and Ea llll can be dis- 
played through the observational parameters of the Seyfert 
1 galaxy 3C 120 and the quasar 3C 345, as shown in Table 

m 

Assuming Qo = 5°, Qo = 60°, 7e = 10, Mo = 1 M© /yr, 
rip ~ 10^" m^'^, j a, I = 0.9, the cross area can be obtained 
as shown in Fig.3. Notice that in Fig.3 the inclined lines 



where 7^ is the Lorenz factor of relativistic electrons, and 
6 = {n — 6)/2, as shown in Fig.2, is the misalignment angle 
between the velocity of protons and the jet axis. 

By Ea ll2l the velocity gained by a baryon at one col- 
lision, is ~ lO^km/s, which is responsible for the velocity 
inferred from the Doppler broadening of emission lines ob- 
served in AGNs. 

Moreover, the jet-disc collision predicts a velocity dis- 
tribution as shown in Fig. 4, in which magnitude of velocity 
decreases with the increase of 6. Such a velocity distribution 
of baryons predicts a bi-polar configuration of line emission. 

In fact, Ea ll2l is obtained ignoring the energy absorp- 
tion by baryons, which can ionize the baryons. For simplicity, 
we can assume that a half of energy before collision is ab- 
sorbed by the baryons, the magnitude of the baryon velocity 
V of Ea ll2l will be smaller. While in order of magnitude the 
conclusion drawn is unaffected. 

Notice that the bi-polar corresponds to half opening 
angle of (tt — 0o)/2, and the smallest defiection angle Oq is 
given by Eq[51 

In fact, Ea ll2l predicts a correlation: the velocity of 
baryons inferred from line broadening of BLR is inversely 
related to precession period of the jet. 
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This is easy to understand. The more the linear mo- 
mentum exchange between electrons of jet and the baryons 
of the disc, the more the momentum and energy gain for 
the baryons, so that they travel faster. On the other hand, 
this corresponds to larger action and reaction forces be- 
tween jet and disc, and hence larger torque on the jet, so 
that the jet precesses more rapidly (and shorter precession 
period). The rate of jet precession can be approximated, 
O ~ Stt/P, which is so small that v of Eg 1121 can be writ- 

1/2 

ten, V ~ y^T^^cQ oc /P. The predicted correlation is 
supported by some AGNs with corresponding observations, 
as shown in Table [2l 



6 DISCUSSION 

Besides interpreting the difficulties in angular momentum 
transfer and precession velocity of jet, a number of observa- 
tions concerning BLR likely favor the new model. 

Firstly, the profiles of the broad lines observed at a high 
spectral resol ution and signal-to-no ise ratio, appear to be 
very smooth (|Eracleous et al.ll2009f ). The velocity distribu- 
tion predicted by our model originates in the jet-disc colli- 
sion, as shown in Fig.4, is automatically smooth. 

Secondly, as shown in Fig.l, the new model actually 
implies a limited jet-disc interaction, in which only a small 
fraction of leptons in the jet is involved in the collision. Once 
the jet plunges into the outer ring heavily, or touches part of 
inner rings, the coUimated jet can be destroyed. This may 
explain the observations of radio quiet quasars, with the 
largest emission line width among AGNs, but weak in radio 
emission. 

Obviously, when a torus is outside the bi-polar config- 
uration of BLR of an AGN , then obscurity of the BLR 
becomes more complicated. Therefore, a LOS with a rela- 
tively small misalignment angle with the jet axis tends to 
observe single peaked emission line of BLR, as some Seyfert 
type 1. 

And if the misalignment between LOS and jet axis is rel- 
atively large, the obscurity depends on the angle of LOS with 
the normal direction of the outer ring. As shown in Fig.l, 
a large misalignment is favorable to observe double-peaked 
BLR, and a small one tends to be obscured, as Seyfert type 
2 (which is absence of BLR emission)^ 

Thirdly, a number of authors (IChen fc Halper"n 1 19891 : 
lEracleous fc Halpernl ll994: Stra tevaet al.ll2006l . booi l have 
pointed out that local viscous dissipation in the line-emitting 
region of accretion disc does not provide enough power to 
account for the observed emission line luminosity. 

Apparently in our model this difficulty is naturally 
avoided. During the collision, as shown in Fig.l, baryons 
at the collision sites (or near) gain energy and momentum 
(linear), which can automatically ionize neutral hydrogen or 
other atoms. Typically, a non-relativistic proton can absorb 
~ leV in a head on collision with an electron of Lorentz 
factor 10 in the jet, which can well account for the energy 
required for line emission of BLR. And multi-collision pro- 
cesses can easily excite a hydrogen to UV band. 

Moreover, if the jet contains certain fraction of baryons 
with a bulk Lorentz factor of 10 also, their collisions with the 
baryons in the outer disc correspond to an energy exchange 



of ~ IKeV, which is responsible for the energy required to 
invoke the line emission of Fe Ka. 

Because NGC 7213 displays very close speeds inferred 
from its BLR and Fe Ka li nes, it is used as an evidence of 
BLR origin for Fe Ka line (|Bianchi et al.ll2008l '). 

From the point of view of the jet-disc interaction, the 
cross section of the relativistic electrons in the jet with the 
non-relativistic Hydrogen in the disc differs from that of 
relativistic protons in jet and non-relativistic iron in the disc. 

Therefore, even if a source shows considerably discrep- 
ancy between speeds inferred from BLR and Fe Ka lines, 
the emission lines may still be originated at the same site of 
collision. So evidence of common origin should be found in 
luminosity, number density, absorbtion extracted from the 
two types of emission lines. 

Consequently, all these issues are realized at the ex- 
pense of a warped disc, which makes the jet-disc interaction 
possible. 
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Figure 1. The interaction of a jet with a warped disc. At the 
collision site, baryons in the outermost ring of the warped disc 
are scattered off, which carry velocity of ~ lO^km/s, and are 
responsible for line emission, as represented by the cloudy con- 
figurations. In the figure, rg is the radius of the outermost ring, 
9f) denotes the opening angle of the jet, oq is the angle of the 
jet with the axis of the precession cone, and ftp is the precession 
angular velocity. Observing at different view angles, the BLR line 
profiles appear different, e.g., 1 is likely responsible for the single- 
peaked one; 3 for double-peak; and 2, 4 for absence of BLR lines 
(obscured by the warped disc or an outer torus co-plane with the 
outer disc). 
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Figure 2. The schematic of collision of a relativistic lepton with 
a non-relativistic baryon. After colliding with a baryon at p. The 
trajectory of the lepton changes for an angle, 0. The exchange of 
momenta between the lepton and the baryon is shown. 
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Figure 3. Parameters space of Vi and ro predicted by the col- 
lision model for sources 3C 120 and 3C 345 respectively. The 
inclined lines correspond to the Vi and ro relation given by Eg llll 
with two different values of M. The vertical ones in the middle 
(solid) correspond to tq given by Eq|3] under two different val- 
ues of M, which correspond to a collision site given by the tidal 
torque. The vertical line at right (dashed) is responsible for the 
radius where the Kepler velocity of gases is equal to that inferred 
from BLR emission lines. The vertical line at loft (doted) cor- 
responds to the disc radius when the Lense-Thirring precession 
yields a precession period of tens of years. 



© 2002 RAS, MNRAS 000, [T]-?? 



The consequence of jet interacting with a warped disc 9 




Figure 4. The velocity distribution of scattered baryons. The 
angle between the velocity of baryons and the jet axis is 6, the 
magnitude of the velocity is given by i; = i;maxCosS, where 
''max = 2y/2^mec/m,p. Note that the velocity field is axial sym- 
metric about the vertical axis, which predicts a bi-polar configu- 
ration of BLR speed. 
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Table 1. Parameters of the Seyfert 1 galaxy 3C 120 and the quasar 3C 345 in the test of jet-disc collision 



Object 


P{yr) 




Ms{Mq) 


Pps (yr) 


rps{pc) 


z 




3C 120 


"12.3 


"3.0 X 10^ 


"'l.O X 10^ 


= 1.4 


=0.002 


0.003 


"1.5° ± 0.3° 


3C 345 


''10.1 


<=4.4 X 10^ 


^=3.6 X 10^ 


^=5.2 


=0.021 


0.5928 


H.3° ±0.5° 



NOTE.-P is the precession period measured in the observer's frame, and Pps is the orbital period of the binary BH system. The mass 
of the primary and secondary BIf are Mp and Ms respectively, rps is the separation between the two BHs, z is the redshift, and ip is 
the half-opening angle of the precession cone given by the model of the following references. 
REFERENCES.-a. Caproni & Abraham 2004b; b. Caproni & Abraham 2004a; c. Caproni & Abraham 2004c. 



Table 2. The correlation of precession period versus the FWHM of Ha emission lines 



Object 


2 


P(yr) 


Pprec (?/^) 


FWHM (km 


Arp 102B 


0.024 


"2.2 


2.148 


''14400 


NGC 1097 


0.004 


"5.5 


5.478 


''7200 


30 345 


0.593 


"10.1 


6.340 


=3600 


OJ 287 


0.306 


"11.6 


8.882 


■^3710 


30 120 


0.033 


"12.3 


11.907 


= 1846 


30 273 


0.158 


"16.0 


13.817 


=3260 


30 279 


0.536 


"22.0 


14.323 


= 1400 



NOTE. -2 is the redshift, and P is the precession period measured in the observer's reference frame P = (1 + 2) Pprec- 
REFERENOE.-a. Oaproni et al. 2004; b. Wang & Zhang 2003; c. Rokaki et al. 2003; d. Decarh et al. 2011. 
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